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Abstract 

We present a new, mathematically rigorous, method suitable for bound state 
and scattering processes calculations for various three atomic or molecular 
systems where the underlying forces are of a hard-core nature. We employed 
this method to calculate the binding energies and the ultra-low energy scat- 
tering phase shifts below as well as above the break-up threshold for the three 
He-atom system. The method is proved to be highly successful and suitable 
for solving the three-body bound state and scattering problem in configuration 
space and thus it paves the way to study various three- atomic systems, and to 
calculate important quantities such as the cross-sections, recombination rates 
etc. 
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The ^He trimer is of interest in various areas of Physical Chemistry and Molecular 
Physics, in particular such as the behavior of atomic clusters under collisions and Bose- 
Einstein condensation. Various theoretical and experimental works have been devoted in 
the past to study its ground state properties and in general the properties of the "^He and 
other noble gas droplets. From the theoretical works we mention here those using Variational 
and Monte Carlo type methods the Faddeev equations and the hyperspherical 

approach P|-pT| . From the experimental works we recall those of Refs. [|T^-|T^ where molec- 
ular clusters consisting of a small number of noble gas atoms were investigated. 

Despite the efforts made to solve the He-trimer problem various questions such as the 
existence of Efimov states and the study of scattering processes at ultra-low energies still have 
not been satisfactorily addressed. In particular for scattering processes there are no works 
which we are aware of apart from a recent study concerning recombination rates [|T6[. There 
are various reasons for this the main one being the fact that the three-body calculations 



*0n leave of absence from the Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, 141980, Russia 



1 



involved are extremely difficult to perform due to the practically hard-core of the inter- 
atomic potential which gives rise to strong numerical inaccuracies that make calculations 
for these molecules cumbersome and unstable. 

In the present work we employed a hard core version of the Boundary Condition Model 
(BCM) |1^,|18[, developed in [jl9|, [^]. The so-called hard-core potentials represent a partic- 



ular variant of the BCM where one requires that the wave function vanishes when particles 
approach each other at a certain distance r = c. Such a requirement is equivalent to an 
introduction of an infinitely strong repulsion between particles at distances r < c. The 
standard formalism for scattering (see, for example, [^) does not deal with hard-core in- 
teractions described by these boundary conditions. Replacement of the finite, for r > 0, but 
often singular at r = 0, repulsive short-range part of the potential with a hard-core inter- 
action turns out to be a very effective way to suppress inaccuracies related to a numerical 
approximation of the Schrodinger operator at short distances. 

In order to outline our method we start from the Schrodinger equation for bound states, 
= E^^ , where \1/ is the three-body bound state wave function. We are concerned with 
states for which E < and that these energies are below the threshold of the continuous 



spectrum of H. Using the Green's formula [^] one can show that the function \l/ satisfies 



the following Lippmann-Schwinger type equation 

^iX) = - J dasGoiX,S;E)^^^{S)-J2 J dX' GoiX, X'; E)iV^^){X') . (1) 



an a 



Here Q is the configuration space of the three-body system in the hard-core model which 
represents only a part of the six- dimensional space, M^, external, |xq| > c^, with respect 
to all three cylinders Fq, Fq = {X G : X = {-Xo,,ya), \^a\ = c^}, a = 1,2,3, where Cq,, 
Ca > 0, stands for the value of |xq,| when the cores of the particles in the pair a contact 
each other. The Xq,, are the usual Jacobi coordinates [^. By Go{X, X'; z) we denote the 
free Green function of the three-body Laplacian —Ax and by n^, the external unit vector 
(directed into Vt) normal to the surface dVL while das is a surface element (five-dimensional 
square) on dVL. By V^, = Kf(xQ,), we denote the pair potentials acting outside the core 
domains, i.e. at |xq| > Cq,. 

The Faddeev components of the function \1/ are introduced via the formulas [|18| , |19[] 

$a(X) = - j dasGo{X,S;E)-^^^{S)- j dX'Go{X,X';E){Vo,^){X'). (2) 
r^rian n 

One can show that they satisfy the following system of differential equations 



.Ax-E)$,(X) = 0, 



(3) 



According to Eqs. (|l|) and @ the sum $q,(X) outside the surface dfl coincides with the total 
wave function ^, i.e., Y^'' <l>p{X) = ^(X), X e . At the same time, it follows from the 



Green's formula that this sum vanishes inside all the core domains, ^1^^ '^isiX) = 0, X e 
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\ f2 . In practice one can replace these, very strong, conditions with the essentially more 
weak ones [|T9| , pO[] 



/3=1 



a 



1,2,3, 



(4) 
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requiring that the sum of $a(X) to be zero only on the cylinders Fq. 

Partial version of the Faddeev equations (0) for a system of three identical bosons read 
(see 



Mil) 



l{l + l) A(A + 1) ^.^ , , 



x^ 



y 



-V{x)'^aL{x,y), x> c 



0, 



X < c 



(5) 



Here, by x, y we denote absolute values of the Jacobi variables x, y and by c, the core size 
which is now the same for all three two-body subsystems. The notation L stands for the total 
angular momentum, and /, A, for the relative angular momenta corresponding respectively 
to a two-body subsystem and complementary particle; a = {/,A}. The potential V{x) is 
supposed to be central, acting in the same way in all the partial waves /. The partial wave 
function \l/aL(x, y) is related to the partial Faddeev components ^aL{x,y) by 



„/ -J -i 



dv haa'{x,y,r]) <^a'L{x',y' 



(6) 



where 



X 



1 2 3 2 

-^x +-y - —xyq. 



' 3 2 ^2 

-X +-y + —xy7], 



with 7] = X ■ y. The explicit form for the function h^^, can be found in Ref. [0,iJ]- The 
functions ^aL{x,y) satisfy the boundary conditions 



^aL{x,y) 1^=0 = ^aL{x,y) = 0, 

while the partial version of the hard-core condition is given by 

'•+1 



duh^^,{c,y,r]) = 0, 



(7) 



(8) 



For the bound-state problem one requires that the functions ^aL{x, y) are square integrable 
in the quadrant x > 0, y > 0. A more detailed and useful in bound state calculations is the 
asymptotic condition 



^aL = ^iJi,uix) hx{^/E - ti^^y) [naL^u + o{y 

V 

exp{i\/Ep + mL/2) 



(9) 



where E is the bound-state energy, p = \X\, 6 = arctany/x, and '4)i^y{x) is the two-body 
partial wave function corresponding to a z/-th bound state e; ^ for the angular momentum 
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value /. The notation h\ is used for the spherical Hankel function. The coefficients aaL,u and 
AaL{(^) describe contributions into from the (2 + 1) and (1 + 1 + 1) channels respectively. 
The corresponding asymptotic boundary conditions for the partial Faddeev components of 
the (2 + 1— >-2 + l; 1 + 1 + 1) scattering wave function as p — > oo and/or ?/ — > oo read as 

+ Y.^v,Ax)h,\^/E^, y) [a5^^;,,(p) + o [y-^')] (10) 



exp(iV£;p + i7rL/2) 



where p = |p| is the relative moment conjugate to the variable y and the scattering energy 
E is given by _E = ei^^+p^. The j\i stands for the spherical Bessel function. The value aj^'^^'^, 
represents, at > eiiy, the partial amplitude of an elastic scattering, a' = a and v' = u, or 
rearrangement, a' ^ a or u' ^ u, process. The functions A^^,'l\6) provide us, at > 0, the 
corresponding partial Faddeev breakup amplitudes. 

We employed the Faddeev equations (||) and the hard-core, (§), and asymptotic, (p|,p!OD, 
boundary conditions to calculate the binding energies of the Helium atomic trimer and the 
ultra-low energy phase shifts of the Helium atom scattered off the Helium diatomic molecule. 
In our calculations we take h^/m = 12.12 KA^. As a ^He-^He interatomic interactions we 
use the HFDHE2 and HFD-B |^6[ potentials of Aziz and co-workers which, we found, 



that they sustain a dimer bound state at — 0.8301 mK and — 1.6854 mK respectively. At 
the same time the ^He atom-^He atom scattering length was found to be 124.7 A for the 
HFDHE2 and 88.6 A for the HFD-B potential. 

In the present work we restrict ourselves to calculations for S-state only. The partial 
components $iAo can be obtained in this case from the addition of even partial waves / and A 
with / = A. The results of the Helium trimer ground-state energy calculations are presented 
in Table I. Although the two potentials used differ only slightly, they produce important 
differences in the ground-state energy. This is in agreement with the finding of Ref. |l^ 



but in disagreement with the statement made in Ref. |Q. It should be further noted that 
most of the contribution to the binding energy stems from the / = A = and / = A = 2 
partial component the latter being more than 35 %. The contribution from the I = X = 4 
partial wave was shown in to be of the order of a few per cent. We have found that the 
Helium trimer can form an excited state with both the HFDHE2 and HFD-B potentials in 
agreement with the findings of Refs. HJ^,!!!!]. Note that in the papers Rllll this state is 



interpreted as an Efimov one ||2^. Our excited state results are given in Table II. 

The phase shift So{E) results, for a Helium atom scattered off a Helium dimer at L = 0, 
are given in Table HI. We considered incident energies below as well as above the breakup 
threshold, i.e., for the (2 + 1 — > 2 + 1) and the (2 + 1 — > 1 + 1 + 1) processes. It is seen 
that, similarly to the bound state results, the inclusion of the / = A = 2 partial wave is 
essential to describe the scattering correctly. 

Our estimation for the scattering length, based on the phase shift results, with the HFD- 
B interactions is 170+5 A in the case where only the / = A = are taken into account and 
145±5 A when both the / = A = and / = A = 2 are considered. We mention here that an 
estimation of igc = 195 A for the ^He atom- ^He dimer scattering length was previously made 
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by Nakaichi-Maeda and Lim P via zero-energy scattering calculations and by employing a 
separable approximation for the HFDHE2 potential. 

The results obtained with two realistic ^He-^He potentials clearly demonstrate the reli- 
ability of our method in three-body bound state and scattering calculations. The effectively 
hard-core inter-atomic potential together with other characteristics of the system, make such 
calculations extremely tedious and numerically unstable. The numerical advantage of our 
approach is already obvious from the structure of Eqs. (|^): When a potential with a strong 
repulsive core is replaced with the hard-core model, one approximates, inside the core do- 
mains, only the Laplacian —Ax, instead of the sum of the Laplacian and a huge repulsive 
term, and in this way a much better numerical approximation can be achieved. Thus the 
present formalism paves the way to study various three-atomic systems, and to calculate 
important quantities such as the cross-sections, recombination rates etc. 
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TABLES 



TABLE I. Bound state energy (in K) results for the Helium trimer. 



Potential 


Faddeev equations 


Variational 
methods 


Adiabatic 
approach 


/ This work H [|7| § 


M M 




HFDHE2 


0.084 0.082 0.092 


0.1173 


0.098 


0,2 0.114 0.107 0.11 


HFD-B 


0.096 0.096 


0.1193 




0,2 0.131 0.130 



TABLE IL Excited state energy (in mK) results for the Helium trimer. 



Potential 


I 


This work 


M 


i 


|11| 


HFDHE2 





1.54 


1.46 


1.04 


1.517 


0,2 


1.74 


1.6 




HFD-B 





2.56 








0,2 


2.83 







TABLE HL The 5-state Helium atom - Helium dimer scattering phase shifts 6q and 6q ' 
obtained with the HFD-B ^He-^He potential. The shifts correspond to the case where only 
the partial wave / = A = was included while the shifts d^^'"^^ were obtained with inclusion of both 
partial waves / = A = and / = A = 2. The values of 6^'"^^ are given in degrees and E, in mK. 



E 


^(0) 

"o 




E 


"o 




E 


"o 


.(0,2) 

"o 


-1.685 


177.5 


177.8 


-1.1 


123.6 


132.6 


0.7 


88.9 


98.7 


-1.68 


172.3 


173.3 


-0.8 


115.0 


124.6 


1.0 


85.7 


95.4 


-1.60 


153.4 


156.4 


-0.4 


105.7 


115.5 


1.4 


81.8 


91.5 


-1.5 


142.3 


148.0 


-0.1 


99.9 


109.7 


1.8 


78.4 


88.0 


-1.4 


135.8 


143.0 


+0.3 


93.9 


103.7 


2.4 


74.2 


83.5 
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